Simple effective interaction: Infinite nuclear matter and finite nuclei 
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The mean field properties and equation of state for asymmetric nuclear matter are studied by 
using a simple effective interaction which has a single finite range Gaussian term. The study of finite 
nuclei with this effective interaction is done by means of constructing a quasilocal energy density 
functional for which the single particle equations take the form of Skryme-Hartree-Fock equations. 
The predictions of binding energies and charge radii of spherical nuclei are found to be compatible 
with the results of standard models as well as experimental data. 

PACS numbers; 21.10.Dr., 21.60.-n., 23.60.+e., 24.10.Jv. 



I. INTRODUCTION 

Last three decades have seen a regular interest to ex- 
plain consistently the properties of nuclear matter, finite 
nuclei and nuclear reactions (nucleon-nucleon, nucleon- 
nucleus and nucleus-nucleus) with an effective interac- 
tion that has the efficiency to describe the two body sys- 
tem accurately. In this context, the study of nuclear 
properties from finite nuclei to highly isospin asymmet- 
ric nuclear matter in a given model is a promising area 
of current interest. Relativistic and non-relativistic mi- 
croscopic models such as Dirac-Brueckner-Hartree-Fock 
(DBHF) 1-9], Brueckner-Hartree-Fock (BHF) [lM3 
and variational calculations using realistic interaction 
pTl [Tsl are considered to be standard references in the 
regime of nuclear matter (NM). The ab initio extension 
of these models to finite nuclei are still in a preliminary 
stage. However, an energy density functional based on 
the microscopic calculations of [16| that reproduces ac- 
curately binding energies and charge radii of finite nuclei 
has been reported recently [l^ . Mean field models using 
effective interactions basically adopt the strategy of fit- 
ting the force parameters simultaneously to finite nuclei 
data and NM constraints (binding energy per particle, 
saturation density, incompressibility, symmetry energy, 
etc.). 

Finite nuclei results provided by Relativistic Mean 
Field (RMF) approaches are quite successful [20l - l23| in 
order to describe the /3-stable and drip-line regions in- 
cluding superheavy elements. Some of these success- 
ful RMF interactions are NL3 [H, NL3* DD-MEl 
[is,], DD-ME2 13 > DD-F [H and DD-MEJ [29i]. In the 
non-relativistic frame the Skyrme [13, HH , Gogny [32l - [36| 
and M3Y [3, [13, [11] effective interactions predict many 
properties of finite nuclei reasonably well. An important 
property of finite nuclei and infinite nuclear systems is 
the neutron and proton effective mass. However, mean 



field models predict neutron-proton (n-p) effective mass 
splittings in isospin asymetric nuclear matter (ANM) not 
always in agreement with the results obtained in micro- 
scopic calculations. For instance, RMF models predict 
that the proton effective mass m* is larger than the neu- 
tron one TO* [39|-[4^, which is contrary to the results of 
microscopic calculations. In the non-relativistic frame, 
the Skyrme (with the exception of the SLy sets [isl - lisi ]) 
as weU as the Gogny forces [32|-l36j. predict n—p effective 
mass splittings which are, in general, consistent with the 
microscopic result m'^>m*. Another important property 
of the mean field models is the mean field generated by 
the effective interaction, which in general is density and 
momentum dependent. In the non-relativistic frame the 
mean fields obtained by means of the Skyrme and finite 
range effective forces differ widely among them. The nu- 
clear mean field corresponding to Skyrme forces has a fc^ 
dependence that is not in agreement with the behaviour 
of the mean field extracted from the analysis of the flow 
data in heavy-ion collision (HIC) experiments at inter- 
mediate and high energies |46l453l | . However, the more 
involved momentum dependence of the mean field pro- 
vided by the finite range forces qualitatively agree with 
the behaviour extracted from the experiment. 

In our earlier works [13, [111 , it has been shown that 
a mean field behaviour consistent with the experimen- 
tal results can be obtained with finite range interactions 
containing a single form factor of Yukawa, Gaussian or 
exponential type. It was also shown [s^ that the range 
of these form factor and the strength of the exchange en- 
ergy in symmetric nuclear matter (SNM) could be con- 
strained using the momentum dependence of the experi- 
mental optical potential extracted from nucleon- nucleus 
collision data [U [1^. In this way the different mean 
fields obtained with the aforementioned form factors have 
a similar behaviour over a wide range of momentum and 
density. This simple effective interaction (SEI) has been 
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used in several studies of NM, such as, momentum depen- 
dence of the mean field and the equation of state (EOS) 
in ANM [HB, IstI . neutron star matter at zero and finite 
temperature 5^, HI] and the thermal evolution of NM 
properties [60 |. 

The main aim of the present work is to show that this 
SEI, able to give an overall good description of NM prop- 
erties in the isospin channel, can also reproduce finite 
nuclei properties with a similar quality to that obtained 
using the more traditional effective interactions such as 
the Skyrme, Gogny and M5F forces as well as with RMF 
parametrizations. Our study of ANM with the SEI, de- 
termines only nine out of the total eleven parameters of 
the interaction. To determine the two remaining param- 
eters is a necessary task for a wider application of the 
interaction. In this work we have adopted the proce- 
dure of fixing them from finite nuclei experimental data. 
Moreover, the knowledge of all the 11 parameters of the 
SEI allows to make additional NM calculations taking ex- 
plicitly into account the spin- and spin — isospin asym- 
metries. 

It is worth mentioning here that full Hartree-Fock (HE) 
calculations in finite nuclei with finite range effective in- 
teractions were performed first using the Gogny force [s^- 
Is^ in a harmonic oscillator basis. Einite nuclei HE cal- 
culations ^ with the 3- Yukawa (M3Y) force [H were 
performed first using the density matrix (DM) expansion 
[62| in the limit of the local Eermi momentum approxima- 
tion [63^ which reduces the HE equations to a local form. 
Later on, both, full HE and Hartree-Eock-Bogoliubov 
(HEB) calculations have been performed with the M3Y- 
interaction [s^, in NM and finite nuclei. An alter- 
native approximation to the HE theory, which is widely 
applied to electron systems, is based in the Kohn-Sham 
(KS) scheme within the framework of the Density 
Eunctional Theory (DET). We are aware that the ap- 
plication of the KS-DET theory to self bound system 
like nuclei is not obvious. However, it has been recently 
shown [sHl that the Hohenberg-Kohn theorem [6^, on 
which the KS-DET scheme is based, has to be reformu- 
lated for the intrinsic frame in the case of nuclei or other 
self quantum systems. The original version of the KS- 
DET theory is local, however, in the nuclear case, non- 
local contributions such as the effective mass and the 
spin-orbit potential are essential ingredients of the en- 



ergy density owing to the momentum dependence of the 
nuclear interaction. The non-local extension of the DET 
was treated first by Gilbert [63|- Later on, a modifica- 
tion of the non-local generalization of the DET and its 
quasilocal reduction were discussed in detail in Ref. [6^ 
and we refer the reader to this reference for more details. 
In order to write explicitly the exchange energy in a local 
form, it is necessary to express the DM in terms of only 
local quantities, such as the particle and kinetic energy 
densities. This can be done, for instance, using the DM 
expansion of Negele and Vautherin [gl] or Campi and 
Bouyssy or, alternatively, with the semiclassical h- 
expansion of the DM [6^. In this work we will use this 
latter approximation together with the SEI to build up 
an energy density functional able to describe ground state 
properties of finite nuclei. It is also important to point 
out that nowadays a modified DM expansion is used in 
Effective Eield Theory calculations with realistic interac- 
tions to write the microscopic energy density in a local 
Skyrme- like form (see (70j and references therein). 

The paper is organized as follows. In Section II we 
present the theoretical formalism for ANM with the SEI 
containing a single Gaussian form factor for describing 
the finite range part. The parameters necessary for a 
complete description of SNM and ANM are discussed 
together with the results. In section HI, a general for- 
mulation of NM having all the three possible asymme- 
tries, namely isospin, spin and spin-isospin, is given in 
order to examine the behaviour of the SEI in the spin 
and spin-isospin channels. In Section IV we build up an 
energy density functional based on the SEI within the 
framework of the quasilocal DET. In the same Section 
the predictive power of the SEI to describe the ground 
state binding energies and charge radii of spherical nuclei 
is analyzed with some detail. Einally, Section IV contains 
a brief summary and conclusions. 



II. ISOSPIN ASYMMETRIC INFINITE 
NUCLEAR MATTER 



Our proposed SEI, that we will use in the calculation 
of NM and finite nuclei properties, has the following form 



Veffir) = to(l + xoP.)Sir) + |(1 + xgP.) (jT^R)) '^^'^^ + + ' ' MP^Pr) fir), (1) 

where, f(r) is the form factor of the finite range interaction which depends on a single parameter a, the range of the 

interaction. In this work /(r) is chosen to be of Gaussian form, 1°' . The other terms have their usual meaning 
54. 15511 . The modified density dependence with the parameter h in the denominator is taken, as in the earlier works 
57H60|, in order to prevent the supraluminous behaviour in NM. The SEI in Eqn. (1) contains 11-parameters, namely 

to, a;o, ^3, a;3, 6, W , B, H, M, 7 and a that will be determined from NM and finite nuclei properties. 

The energy density functional HxiPm Pp) in ANM at temperature T, derived using the SEI in Eqn. ([l}, can be 

expressed in terms of the neutron (proton) densities p„ (pp) and their respective momentum distribution functions 



3 



/^(k) {fP{k)) as M: 



2m J 



1 
2 

_^ 



fo 



„7+l 



l + bp 



[pI 



£o_ 
Po 



„7+l 

Po 



P 



l + bp 



PnPp 



2pi 



where, po is the saturation density of SNM and gex(k.) 
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is the normahzed Fourier transform of the finite 



range form factor f{r). The index I [ul) is used for the interaction between a pair of hke (unhke) nucleons. For the 
sake of simphcity, we have assumed that the interaction between a hke and unhke pair of nucleons have the same 
range but differ in strengths. The energy density in Eqn. ^ is written in terms of nine parameters, namely, 7, b, 
Eq, Eg', e5y,e"', eia;, £"2: a. The six new parameters are connected to the interaction parameters of Eq.([T]) through 
the relations 
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The neutron (proton) single particle energy e^^^''(k, p„, pp) can be obtained as the functional derivative of the energy 
density _ffy(p„, pp) in Eqn. ([5]) with respect to the neutron (proton) and is given by 
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where, the first term corresponds to the kinetic energy of neutrons (protons) and u^*^^-* is the corresponding mean 
fields. 

The Fermi-Dirac momentum distribution functions fxO^) ^-nd /-^(k) at finite T take the form of step- functions in 
the T=0 limit. In this case the mean fields u"*^^'^ (k, p„, pp) and energy density -ff (p„, Pp) can be evaluated analytically. 
The corresponding expressions for a Gaussian form factor /(r) read, 
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Here p = p„ + is the total nucleonic density, fc„(p) = (37r^p„(p)) ^ is the neutron (proton) Fermi momentum and 
A = 2/a. The exphcit expression for the energy density and other relevant quantities in ANM obtained with a Yukawa 
form factor can be found in Refs. [57l - l60| . It is worth mentioning here that compact formulae for energy in NM for 
several finite range form factors can be found in Appendix A of Ref. [s^] ■ 

In the limit of symmetric nuclear matter, Pn = Pp = p/'2, the corresponding expressions H{p) and u(k, p) become. 



H{p) = pe{p) 



SfPkjp 



10m 

\^ex ' ^ ex J 



4/90 



-P 



4po 
3A6 9A4 



4p2 



,7+1 



p(R) 



1 + V(R) 



16A:6 8fc4 



3A4 



3A6 \ 



^ 2kl 



3^3 r2k,JK 



dt 



(7) 



and 



u(k,p) 



2/90 



2po 



-P 



2pr 

3A4 



' k^k £ \ 2 



3A3 

4fc| J(^z^) 



(8) 



/ 2 

where fc/ = l^f-p) is the Fermi momentum and e{p) is the energy per particle in SNM. The new parameter 



combinations, 



(9) 



alongwith 7, 6 and a are the six parameters needed for 
a complete description of the mean field properties and 
the EOS in SNM. The way of determining the parame- 
ters which enter in Eqs. Q and ^ in this case is similar 
to the one described in Ref. [58| for a Yukawa form fac- 
tor. The only difference in this work is that we take the 
kinetic energy term in its non-relativistic form for conve- 
nience of application to finite nuclei. The range parame- 
ter a and the exchange strength combination Sex in SNM 
are determined by adopting a simultaneous optimization 
procedure with the constraint that the attractive opti- 
cal potential changes sign for a kinetic energy ^|^=300 



MeV of the incident nucleon (see Ref. [13] for details). 
In this way the values eea;=-94.46 MeV and a=0.7596 
fm are found where we have used only the standard 
values of the nucleon mass m=939 MeV, Fermi kinetic 
energy h'^k'j^/2m=36A MeV (the Fermi momentum is 

kfg = {3n'^po/2y^^ with a saturation density po=0-157 
fm~^) and energy per particle in SNM e{po)= -16.0 
MeV. The momentum dependence of the nuclear mean 
field in SNM computed using this SEI for three different 
densities, p=0.1, 0.3, and 0.5 fm~^, is displayed in the 
upper panel of the Fig. 1. In order to compare with the 
predictions of the realistic interaction UV14-I-UVII [l8| . 
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we plot u'^^Ck^p) = [^(k, p) — u(kf,p)] as functions of 
k. Here the saturation properties of the interaction have 
been substracted because they may be model dependent. 
We find a quite good agreement between the results com- 
puted with the SEI and the microscopic predictions over 
a wide range of momenta and densities. The same com- 
parison is shown in the lower panel of Fig 1. for Gogny 
Dl ^2], DIS 34] and DIM 36] sets with the microscopic 
results 18]. We see that the momentum dependence of 
the Gogny Dl force agrees well with the microscopic re- 
sults upto density p = 0.3 fm~^ over the whole range 
of momentum shown in the figure, but agreement of the 
results become qualitative at higher density as can be 
seen from the comparison for the curve corresponding to 
p = 0.5 fm~^ . On the otherhand, the Gogny DIS and 
DIM forces do not follow the trends of the microscopic 
calculation beyond k i 2 fm~^. It may be pointed out 
that the Gogny forces are fitted basically to finite nuclei 
data and, therefore, there is no reason a priori for re- 
producing the momentum dependence of the microscopic 
calculation in NM over a wide range of momentum in case 
of all the Gogny force sets. 

The effective mass in SNM is given by. 
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which is momentum and density dependent and can be 
calculated once the exchange part of the mean field in 
Eqn. ([5]) is known. For the values of the exchange 
strength and range determined above, the prediction of 
the effective mass at saturation density p — po and mo- 
mentum k = fc/p is m* /m=0.709. The parameter b has 
been adjusted to avoid a supraluminous behaviour at 
high densities in SNM at T=0 (see Ref. [7l| for more 
details). With the standard values mentioned before we 
obtain 6=0.5914 fm^. The two remaining strength pa- 
rameters eo and are obtained from the saturation con- 
ditions. The exponent 7, that determines the stiffness of 
the EOS of SNM, is taken as 7=| which gives a value 
of incompressibility at normal NM density, K{po)=245 
MeV. With the six parameters of SNM determined, one 
can calculate the energy per particle e{p) and pressure 
P{p) in SNM. These quantities are displayed as a func- 
tions of the density in Fig. 2(a) and 2 (b), respectively. 
They are compared with the values of e{p) and P{p) pre- 
dicted by some microscopic calculations and several non- 
relativistic and relativistic models. In Fig. 2(b) we also 
display the band of allowed values of P{p) in the range 
2 — 4.6 fm~^ extracted from the analysis of the flow data 
in high energy HIC experiments [72|- Further experi- 
mental information of the pressure-density relationship 
is provided by the analysis of production data [73] 
in the low density domain between 1.2 — 2po fm~^. We 
find that our SEI curve P{p), obtained using a value of 
7=1/2 passes well within these experimentally extracted 
regions. 

The complete study of ANM now requires the split- 
ting of the three strength parameters e^x = "° 9 "° i 



£0 — ^ and £^ = '^"'~!^^"' into two specific channels 

for interactions between like and unlike pairs of nucleons. 
We have fixed the splitting of £ex into the like channel 
to be Sgj. = 2eea:/3. This is the critical value for which 
the thermal evolution of NM properties in pure neutron 
matter (PNM) does not surpass the SNM results at any 
density for any temperature [13] • The splitting of the 
exchange strength parameter Eex into {I) and {ul) chan- 
nels fully determines the effective mass behaviour in 
ANM. The effective masses in ANM are calculated 
starting from their usual definitions. 
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The magnitude of the splitting of the n,p-effective masses 
in ANM at normal density, [m*/m]„ — [m*/m]p, obtained 
with the SEI is shown in Fig. 3 as a function of the 
isospin asymmetry /? = ^^"^'^p-' alongwith the results 
of microscopic calculations and the values obtained with 
the Gogny effective forces. Our result compares well 
with the prediction of the DBHF calculations using the 
Bonn potential [3] over the whole range of asymmetry /3. 
The Gogny Dl interaction also predict closely the sim- 
ilar trend as that of our SEI. But, the n — p effective 
mass splittings obtained using the Gogny DIS and DIM 
interactions are very similar between them and have a 
smaller value as compared with DBHF result. On the 
other hand, the splittings in the cases of BHF-I-3-BF and 
EBHF-I-3-BF [1, [1, [74] calculations are found to be close 
to each other, predicting a larger n — p effective mass 
splitting than the DBHF model. 

Once the like and unlike components of Sgx are fixed, 
the splitting of the two remaining strength parameters 

combinations, namely eo = and e-y ~ ("^t+^t ) 

can be obtained from values of symmetry energy Es{pq) 
and its derivatives E'^{po) = Po^^^i^ computed at satu- 
ration density. We have assumed a value of Es{po)—35.0 
MeV which is within the possible range between 30-35 
MeV [3l|,[7i-[7i|. The value of S^(po) is ascertained from 
an universal high density behaviour of the asymmetric 
contribution of the nucleonic part of the energy density in 
charge neutral beta-stahle n+p+e+v matter |28l . [58l |. The 
value E'^{pq) =25.42 MeV obtained in this way predicts 
a slope parameter (symmetry pressure) L(po)=76.3 MeV 
that is well within the range 70 ± 15 MeV of recent Fi- 
nite Range Droplet Model (FRDM) prediction Our 
value L{pq) also lie within the window 45-80 MeV esti- 
mated from a compilation of different predictions of L 
deduced from experiments involving antiprotonic atoms, 
heavy-ion reactions, proton scattering, nuclear masses, 
microscopic calculations and giant dipole resonances [s^] . 
The values of the nine parameters altogether necessary 
for a complete description of SNM and ANM alongwith 
the NM properties (saturation density po, binding en- 
ergy per nucleon in SNM e{po), incompressibility K{po), 
symmetry energy Es{pa) and slope parameter L{pq)) are 
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given in Table I. The strength parameters in SNM can 
be obtained from this table using the relations of Eqn. 

Once the nine parameters that fully determine the 
ANM are known, the energy per particle in PNM, (p), 
can be easily calculated and their values as a function of 
the density p are displayed in Fig. 4 alongwith the results 
of realistic 81] and some effective models. In the low den- 
sity region, p < 0.12po, it has been verified that our curve 
passes well through the region predicted by different mi- 
croscopic calculations performed using Quantum Monte 
Carlo techniques [82| as it can be seen in the inset of this 
figure. In the high density region the SEI curve shows 
a stiff behaviour because of the relatively high value of 
Es{p). The effective M3Y-P5 [ss] curve has a relatively 
more soft behaviour throughout the density range but 
maintains an increasing trend unlike the case with the 
DIS Gogny force displayed in the same figure. Further, 
the EOS of PNM predicts a ratio 3P^(po)/(ipo) =1-0 
MeV /m^"^, which is in agreement with the value ob- 
tained by Piekarewicz [s^] from an analysis of the rela- 
tion between the slope parameter L and the pressure of 
PNM at normal density po. The density dependence of 
nuclear symmetry energy, Es{p), calculated as the differ- 
ence between the energy per particle in PNM and SNM, 
is displayed in Fig. 5. It compares well with the results 
provided by the realistic Al8+(5v-|-UIX* interaction flT\ 
upto a density 0.2 fm^^ and follows a similar trend till 
a density 0.8 fm^^. Our curve also compares reasonably 
well with the DBHF [1, Q results upto a density of 0.5 
fm~^, but clearly differs beyond this point where this 
microscopic calculation shows a very stiff behaviour of 
the symmetry energy with increasing density. 

III. NUCLEAR MATTER UNDER ISOSPIN, 
SPIN AND SPIN-ISOSPIN ASYMMETRIES 

In the previous section we have restricted our NM 
study to the case of isospin asymmetric NM, i.e., ANM, 
which has much physical relevance from the point of view 
of the available experimental and empirical constraints. 
However, in order to examine the predictions of SEI in 
the spin and spin-isospin channels, we shall make first the 
explicit formulation of NM having both spin and isospin 
asymmetries. The total NM density is now given by 

P ^ Pnu + Pnd + Ppu + Ppd, (12) 

where, nu, nd {pu, pd) denote the neutron spin-up and 
spin-down (proton spin-up and spin-down) states and 
Pnu, Pnd, Ppu and ppd are the respective densities. The 
Fermi-momenta corresponding to these four densities are 
expressed as 

= (e^V^y/', (13) 

with i=n,p and j=u,d. Making Taylor expansion of 
the energy density around the spin-saturated SNM value 



H(p) and keeping terms upto lowest order only, one can 
obtain 

H{pnu,Pnd,Ppu,Ppd) = H{p) + ^alpErip) 

+ ^alpE,{p) + ^al^pE„r{p), (14) 

where ar, and a^r are the isospin, spin and spin — 
isospin asymmetries defined as 
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(15) 



and £'r, E^ and E^r are the respective symmetry en- 
ergies. The expressions of these symmetry energies are 
given by, 

(16) 

where, ji, I — 0,1 are the spherical Bessel function of or- 
der I, Xf = akf, kf being the Fermi-momentum in spin 
saturated SNM, and a is the range of the interaction. In 
these equations w"^, v*'^ , and w"" are the interactions 
in the singlet — even, triplet — even, triplet — odd and 
singlet — odd states respectively. Calculation of the sym- 
metry energies using Eqn. (|16|) requires to know the in- 
teractions in all these four different states independently. 
In the case of the SEI, this is possible provided that all 
the eleven parameters of the interaction are known. The 
study of ANM performed in the previous section allows 
to fix nine of them. The two open parameters {to and 
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xq considered here) are determined from the study of fi- 
nite nuclei as wih be explained in the next Section. To 
complet the study of NM matter in the spin and spin- 
isospin channels we will use the numerical value of the 
parameters of the SEI interaction reported in Table II. 

We first calculate the contributions of the singlet — 
odd (SO), singlet — even (SE), triplet — odd (TO) and 
triplet — even (TE) states to the potential energy per 
nucleon {V)/A in spin saturated SNM. These contribu- 
tions, computed using the SEI, are displayed in the four 
panels of Fig. 6 as functions of the Fermi-momentum kf 
along with the results for the Gogny DIS and M3Y-P5 
interactions. The odd-states contributions of SEI and 
DIS have a similar behaviour, being attractive in the SO 
and repulsive in the TO channels. However, these contri- 
butions computed with the M3Y-P5 force show a rather 
strong repulsive character in both odd channels. In the 
TE channel the contributions computed with the three 
interactions have a similar behaviour, being attractive 
upto kf ^ 2 fm^^ and thereafter they become strongly 
repulsive. In the SE channel, the contributions computed 
with the three forces show the same behaviour upto a 
Fermi-momentum corresponding to normal NM density. 
Thereafter the attraction decreases and becomes repul- 
sive at relatively small kf value in case of SEI and at 
a larger kf value for M3Y-P5. On the contrary, the at- 
tractive contribution of DIS in the SE channel increases 
when the Fermi momentum kf increases. 

The density dependence of isospin, spin and spin — 
isospin symmetry energies computed with Eqns. p6p us- 
ing the SEI are displayed in the three panels of Fig. 7. 
The corresponding results obtained with the Gogny DIS 
dl], DIN [3i| and DIM ^ forces are also shown in the 
respective panels. The isospin symmetry energy Er, in 
case of the SEI, has a strong repulsive density dependence 
in comparison with the considered Gogny forces and does 
not predict a preferred PNM than SNM at any density. 
In the case of the spin symmetry energy E^, however, 
the Gogny forces show a very strong density dependence 
as compared to the behaviour exhibit by the SEI. This 
is mainly due to the 2:3 = 1 value in the Gogny forces, 
which makes the contribution of the density dependent 
term of the spin symmetry energy highly repulsive. In 
the spin channel although the SEI has a relatively soft 
behaviour, as the DIM Gogny force in isospin channel, 
it does not exhibit instability related to the transition of 
unpolarized to polarized PNM. A value of Edpo)— 26-55 
MeV is obtained in comparison to 27.57 MeV (Dl), 29.13 
MeV (DIS), 22.69 MeV (DIN), 28.73 MeV (DIM) and 
41.01 MeV (M3Y-P5). The spin — isospin symmetry 
energy Era- shows a similar trend either computed with 
the SEI or with the Gogny forces, except Dl. It reaches 
a maximum value and decreases becoming zero at cer- 
tain densities and thereafter remains negative. However, 
the Gogny Dl interaction shows extra stability in the 
spin — isospin channel. The density at which the curve of 
SEI crosses the x-axis is larger than the crossing density 
in the case of the three curves corresponding to Gogny 
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k (fm"^) 

FIG. 1: (Color Online) The momentum dependent part of 
the mean field u^^{k, p) as a function of momentum k for the 
SEI at three different densities, p=0.1, 0.3, and 0.5 fm~^, 
are compared with the predictions of the realistic interaction 
UV14-I-UVII [13] in the upper panel and the same comparison 
for the Gogny Dl [H, DIS [11] and DIM Q forces in lower 
panel. 



DIS, DIN and DIM. A value of £;^^(po)=23.09 MeV is 
obtained as compared to 30.13 MeV (Dl), 29.13 MeV 
(DIS), 22.69 MeV (DIN), 28.73 MeV (DIM) and 41.01 
MeV (M3Y-P5). In Fig. 7(a), the symmetry energy, 
calculated as the difference of the energy per particle in 
PNM and SNM, Es{p) = [e^(p) - e{p)] is also shown 
for comparison with the corresponding results of E^ip) 
obtained from the Taylor series expansion of the energy 
density in Eqn. p^ . In NM- calculations, the former def- 
inition of the symmetric energy is widely used. It is exact 
at the two extremes of isospin asymmetry. The compar- 
ison between the two, Es and Er, is excellent over the 
whole range of density. The small difference is attributed 
to the higher order contributions of the Taylor series ex- 
pansion which have been neglected. 



IV. FINITE NUCLEI 

As we mentioned at the Introduction, our aim in this 
Section is to explore the ability of the SEI in describing 
ground-state properties of finite nuclei. In Section II we 
have shown that information of NM allows to determine 
nine of the eleven parameters of the SEI. Keeping these 
nine parameters fixed, we will determine the two remain- 
ing open parameters, namely to and xq, as well as the 
strength of the spin-orbit interaction, Wq, from experi- 
mental data of some magic nuclei. Once the SEI is fully 
determined, we will perform some additional calculations 
to check the predictions of our model in other finite nu- 
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TABLE I: The nine parameters of ANM for the simple effective interaction (SEI). The connection of the new parameters with 
the parameters of SEI is in Eqn. (3) . 
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b 
fm 


a 

fm 


MeV 


ul 

MeV 


MeV 


MeV 


4 

MeV 


£o 
MeV 


1 

2 


0.5914 


0.7596 


-94.76 


-125.95 


77.507 


97.247 


-78.78 


-111.69 


Nuclear matter properties at saturation condition 


Po ifm 


e{po) (MeV) 


K(po) (MeV) 


l^(PO,fc/o) 


Es(po) (MeV) 


L{po) (MeV) 


0.157 


-16.0 


245 


0.709 


35.0 


76.26 




FIG. 2: (Color Onhne)(a) Energy per particle in SNM as a 
function of density for the SEI compared with DIS ^34), M3Y- 
P5 i3a], DBHF \g, NL3i2a], G2 ^] and realistic calculation 
[itI ]. (b) Pressure as a function of density in SNM for SEI 
compared with DBHF g, i], DD-F [13 , HIC [zl and k+ 
production data [t^]. See text for details. 
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FIG. 4: (Color Online) The neutron matter EOS for SEI com- 
pared with DIS [13 , M3Y-P5 [38| and Realistic interaction 
[Slj . The behaviour in the low density region, p/po < 0.12, is 
shown in the inset figure, where the Monte Carlo simulation 
results and results of microscopic calculations are compared 
(see text for details). 
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FIG. 3: (Color Online) The neutron and proton effective 
mass difference in ANM at normal density po as a function 
of isospin asymmetry /3 for SEI compared with DBHF 
BHF-H3-BF, EBHF+3-BF [1,1, [H and Gogny Dl [H, DIS 
[13 and DIM |3^ sets. 
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FIG. 5: (Color Online) Density dependence of the sym- 
metry energy for SEI compared with realistic interaction 
A18+5t)+UIX* [i3l and DBHF gi]. 
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FIG. 6: (Color Online) Contribution of the SO, TO, SE, and 
TE channels to potential energy per paricle in spin saturated 
SNM. The predictions of Gogny (DIS) [si] and M3Y-P5 ^ 
are also shown for comparison. 




P (fm" ) 



FIG. 7: (Color Online) The isospin, spin and spin — isospin 
symmetry energies are shown as functions of density. The 
corresponding results for Gogny (Dl) [H, (DIS) H, (DIN) 
[ssj and (DIM) [s^ sets are also shown in the respective fig- 
ures. In (a) the result of Es(p) denoted by SEI(_Es) is shown 
for comparison. 



clei. We will restrict this preliminary study to spherical 
nuclei only and leave for a future work the investigation 
of the deformation properties of the SEI. For open shell 
nuclei we have to add pairing correlations. As we will ex- 
plain below with some detail, we use the BCS approach 
together with a zero-range density-dependent pairing in- 
teraction that was devised to simulate the behaviour of 
the Gogny Dl pairing gap in neutron matter [sij . 



A. The quasilocal Density Functional Theory 

To deal with finite nuclei, we apply the KS 16411 method 
within the framework of the non-local DFT 63. As it 
is proved in this reference, the Lieb theorem 85|, which 
establish the many-to-one mapping of a A-particle Slater 
determinant wavefunctions onto the local particle 
density p(r), allows to write the energy density func- 
tional in the non-local case as e[po] — eo[po] + Erc[p]- 
Here po is the Slater determinant DM with occupation 
numbers either 1 or and Sq has the form of the HF 
energy obtained with an effective Hamiltonian H. The 
remaining part, Encip], that is a functional of the lo- 
cal density only, is the residual correlation energy. This 
contribution accounts for the difference between the ex- 
act energy functional, constructed with the true micro- 
scopic Hamiltonian H, and the HF energy provided by 
the reference Hamiltonian H. In Ref. [g^ it is also shown 
that another reduction can be performed by mapping 
the Slater DM onto a set pql of local particle, kinetic 
energy and spin densities for neutron and proton, i.e. 
pQL ^ (/9„, pp, r„, Tp, Jn, Jp). The particle, kinetic en- 
ergy and spin local densities entering in the set p'^^ are 
obtained from the single-particle orbitals that define 
the Slater determinant 'So as 



(17) 



(18) 



and 



J9(r) =i^^0*(r,tT,g) [{a)a,a' x V] </.,(r, a, g). (19) 



respectively. 

Using this reduction, it is possible to write finally the 
energy density functional in a quasilocal form as e[p'^'^] = 
£o[p'^^] + Ejic[p]- By applying the variational principle 
to the functional e[p^^] and using as functional variables 
the single-particle orbitals (f> and 4>* , one obtains the cor- 
responding KS equations. We shall point out that within 
this quasilocal theory, one would be able to compute the 
exact ground-state energy and the exact local particle 
densities if the exact density functional elp'^^] is known. 
However, the kinetic energy densities, t„ and Tp, as well 
as the spin densities, Jn and Jp, correspond to uncorre- 
lated system and therefore do not concide with the exact 
densities within this approach. As explained in [68j . there 
is some freedom in the choice of the effective Hamiltonian 
H. Therefore, we can choose H as an ^-particle effective 
interaction of the type H = T + J2i^j vf^^ + Y.i^j 
The nucleon-nucleon interaction wf^^ contains the con- 
triutions of the effective force that we take in this work 
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as the density-independent finite-range part of tlie SEI 
in Eqn. ([ij, and tlie spin-orbit contribution cfiosen in 
the form used in the Skyrme and Gogny interactions: 
vff = iWo (cTi -I- cTj) [k' X (5(ri,rj)k]. From this effective 

Hamiltonian H we obtain the quasilocal energy func- 
tional as: 

[p'^^] = juod^R, (20) 
where the energy density T-Lq reads 

/j2 ^^^^ ^^^^ 

(21) 

The kinetic energy part corresponds to the non- 
I 



interacting contribution given by Eqn. ([18]). The 
Coulomb energy is taken in the usual way as the direct 
term plus the exchange contribution computed at Slater 
level using the point proton density: 

1 

The spin-orbit energy density, computed using the afore- 
mentioned zero-range force, becomes: 

H^O(R) = _ l^^g [p(R)VJ + p„(R)VJ„ + pp(R)VJp](23) 

The direct contribution to the nuclear energy H^""^' com- 
ing from finite range part of the SEI is given by 



Nucl 



)p(ri)p(r2) 



M 

-H" + Y ) [Pn(ri)Pp(r2) + Pp(ri)p„(r2)] /(|ri - r2|) 



(24) 



All these local contributions to the energy density Ho constitute the so-called Hartree part of the functional. Up 
to this point, we have developed the exact theory. In the next step we shall make some approximations, similar to 
those used in Refs. [Hll^lS^. To compute the quasilocal energy density corresponding to the exchange terms of the 
nucleon-nucleon force, SEI in this case, we use the Extended Thomas-Fermi (ETF) expansion of the DM up to 
order which is widely discussed in Ref. [63| • For spin-saturated nuclei the exchange nuclear energy density splits into 
two part 



rL. 



Nucl 
exch 



n 



Nucl 
exch 



+ "^, 



Nucl 
exch,2' 



(25) 



The first term corresponds to the zeroth order of the fi-expansion (Slater approximation for the density matrix) given 

by 



n 



Nucl 
exch.O 



2 2 



3ji(fcgr) 



M - 



H 



/r,N3jl(fcn7;) ,„.3ji(fcpr) 



(krir) 



(kpr) 



(26) 



where, r = ri — r2 and R = '"^^'"^ are the relative and center of mass co-ordinates, respectively. In Eq. ([26l) fcg(R) = 
[Sn^PqCR)]^/^ and ji{x) is the _^herical Bessel function. The second term of Eqn. ([25]) is the fi^— contribution to the 
exchange energy which reads josl l69l | 



n 



Nucl 
exch,2 



E 

q=n,p 



2m 



i9q - 1) 



-y^Pq 



(27) 



Here gq = ^^(R, fcg) and g'^ = {dgq{R,kq)/dk)k=kg. 
momentum dependent effective mass given by 

g,(R,fc) = l 

In this equation Vg^"^'^ is the Wigner transform of the exchange potential. 



V27 pq 36 
The function gq — gq{'R,k) is the inverse of the position and 



dk 



^exch,q(^^ k) 



H 



W 



M+--B~- — ]pqiR 



(kqr) 



M+^\pq,{R) 



.^.hikq'r) 



{kq,r) 



(28) 



(29) 



where q = n,p and q' — p, n. Notice that in the ANM limit Vj^^^^^CR, k) is just the exchange part of the mean field 
M'(k, Pq, pqi) in Eqn. ^ and therefore, in this limit, the inverse effective mass gq — (7g(R, kq) reduces to Eqn. (jlip . 
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It is worth noting that within the semiclassical ETF expansion of the DM the kinetic energy is a functional of 
the local density only. However, it was found in |69l] that the use of the quanta! kinetic energy density (ITS)) in the 
/i^-contribution to the exchange energy (f27| . significantly improves the agreement with the corrersponding full HF 
calculation. Therefore, we will use this ansatz here as it was done in previous works [gs', '69', '89'] . It is also important 
to note that we have repalced the exact quasilocal functional Eq [p*^^] by an approximated one calculated with the 
ETF prescription. The difference between them gives a very small contribution that cannot be completely included in 
the residual correlation energy because the difference depends on p^^ (through r„, Tp, J„ and Jp) while -E_rc depends 
on p only. Notice, however, that this small difference, due to the localization of the exchange energy, persists if one 
uses another DM expansions as the Negele-Vautherin [g^ or Campi-Bouyssy [63| ones. 

Let us now discuss the residual correlation energy Ej^c part of the total energy density. The SEI, as it happens with 
another finite-range forces as the Gogny or M3Y ones, can be split into a density-independent and a density-dependent 
parts. Therefore, a reasonable ansatz, in the spirit of the DFT, is to take the residual correlation energy Ej^c as the 
HF contribution provided by the density-dependent part of the interaction in ANM in a local density approximation: 

Erc = [(1 - 2;o) [pUR) + plm] + (4 + 2xo)p„(R)pp(R)] d^R 

+ 1 / [(l-a;3)[p^(R)+p^(R)] +(4 + 2:.3)p„(R)pp(R)] (^-^^^^J d^R. (30) 

As mentioned before, the variational principle applied to the full functional e[p'^^] = EqIp'^^] + Eiic[p] allows to 
obtain the following set of KS single particle equations: 

hq4>i = £i(j)z, (31) 

where, 

/ig = -V^^J^V + C/,(R)-»W,(R).[Vxa], (32) 

and 

f,„(R,.#:i. W,,R,.#^. (33, 



2m* (R) 9r,(R)' ' dpq{R)' ' 9Jg(R)' 



that are formally similar to the equations of motion ob- 
tained with zero-range Skyrme forces. The quasilocal 
DFT presented in this work is, actually, very similar to 
the one used by Hoffman and Lenske in Ref. |2j]. The 
main difference is that the localization of the exchange 
energy is performed here using the ETF expansion of the 
DM while in Ref. 2] the Campi-Bouyssy expansion of 
the DM is used. As it is pointed out in this reference 0, 
the DM expansions retain, to certain extension, the non- 
local effects of the underlying interaction. In the case 
of the ETF expansion of the DM these effects are col- 
lected in the effective mass and its derivative respect to 
the momentum which appear in Eqn. [571 

In Ref. [6^ we have chosen the parameters of wjj^ in H 
and the ones of Ejic to be equal to the parameters of the 
Gogny DIS force. In this way we check to which extend 
the quasilocal KS-DFT approach is able to reproduce the 
full HF results with this interaction. We have performed 
the same comparison with the DIN and DIM forces. It is 
found that the difference between the HF and the quasilo- 
cal KS-DFT energy per nucleon in magic nuclei is always 
smaller than 0.01 MeV in ^^Ca, ^^Zr and ^"Spb for all the 



Gogny forces considered. For light nuclei the differences 
in the binding energy per nucleon are a little bit larger 
but smaller than 0.09 MeV per nucleon in ^^O and than 
0.05 MeV per nucleon in "'"Ca for these Gogny forces (See 
also Table I of ) . These results show that the quasilo- 
cal energy density functional obtained with the KS-DFT 
formalism explained here describes ground-state proper- 
ties of finite nuclei fairly well. This fact motivates us 
to use directly the quasilocal energy density functional 
e[p'^-^] obtained with the SEI fitting its undetermined 
parameters, to and xq of the E^c and the strength of 
the spin-orbit interaction Wq, to experiemntal data. 

Certainly, pairing correlations play an important role 
for open shell nuclei. To deal with such nuclei, it is 
mandatory to include the pairing correlations. The first 
formal generalization of the Hohenberg-Kohn theorem to 
paired systems was performed in Ref. [86] for supercon- 
ductors. Some modifications of this approach were devel- 
oped later in Refs. [13, H^]. More recently, in Ref. [11], an 
extension of the DFT including pairing correlations with- 
out formal violation of the particle-number conservation 
and the quasilocal reduction of this non-local theory are 



12 



TABLE II: The parameters of the simple effective interaction (SEI) and spin-orbit strength Wo- 
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MeV 
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0.5914 


437.0 


0.6 


9955.2 


-0.118 


-589.09 


130.36 


-272.42 


-192.16 


0.7596 


115.0 



TABLE III: The calculated binding energy per particle BE/A, charge radius r^h, the root mean square radii for neutr on r„ 
and proton are compared with the Gogny (DIS) [Hi, M3Y-P5 [11], RMF (NL3*) ^ and experimental data [lo3,[io2|. The 
energy is in MeV and radius in fm. 
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7.867 
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Gogny (DIS) 
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M3Y-P5 
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discussed with detail. The equations of motion associated 
to this energy density functional including pairing corre- 
lations have the same form as the HFB equations. How- 
ever, for pairing calculations of nuclei not too far from 
the /3-stability line, the simpler BCS approach can be suf- 
ficient for describing their ground-state energies [90l | . For 
these nuclei, the Fermi level lies appreciably below zero 
and, consequently, the levels around it, which mainly con- 
tribute to pairing correlations, are also well bound avoid- 



ing the problems of the standard BCS approach near the 
drip lines 0, [o^]. For practical BCS calculations, one 
extracts from the energy density functional a part which 
depends only on the normal density and that describes 
the nuclei without pairing effects. The remaining part 
contains the contributions arising from the pairing cor- 
relations. In the BCS approach the motion equations, 
in the case of spherical symmetry, reduces to a set of 
mean-field single-particle equations, that in the quasilo- 
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FIG. 8: The calculated single particle energy levels for 
are compared with the Gogny (PIS) HI, M3Y-P5 [H 
the experimental data [l04l . llOSi ]. 



and 



TABLE IV: The rma of Sin and Sn deviations for SEI com- 
pared with the corresponding results obtained in other HFB 
calculations. 



Z or N SinISn 
Chain SEI 



S2N /Sn 
Ref. fllll 



S2n/Sn 

Ref. [1121 



Z= 20 0.85/0.62 
Z= 50 0.55/0.62 
Z= 82 0.46/0.82 
N= 50 0.22/0.24 
N= 82 0.21/0.42 
N= 126 0.67/0.51 



0.82/0.76 
0.29/0.21 
0.23/0.37 
0.37/0.26 
0.43/0.31 
0.42/0.23 



1.02/0.92 
0.43/0.35 
0.58/0.53 
0.41/0.23 
0.50/0.56 
0.88/0.52 



cal approach are given by Eq. (|31l) . and another set of 
gap equations 

47r 



A,; 



4Tr 2Ek' 



E, 



(34) 



In the gap equations (Eqn. (1341) ). arc the quasipar- 
ticle energies with £i the single-particle energies, the 
state dependent gaps and A the chemical potential to en- 
sure the right number of particles. In this equation V^jf' 
are the reduced matrix elements of the effective inter- 
action in the particle-particle channel. Eqns. (I3ip and 
(l34t are coupled among them and have to be solved self- 
consistenly. In Refs. [SS, |93| the ability of this quasilocal 
extension of the DFT including pairing correlations at 
BCS level to reproduce full HFB binding energies and 
gaps was checked. The calculations were performed for 
several Sn and Pb isotopes lying in the stability valley us- 
ing the Gogny Dl and DlS forces in both, particle-hole 
and particle-particle, channels. In these test calculations 



TABLE V: The rms of S2N deviations for NL3, SLy4, SkM* 
and DlS 
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FIG. 9: The density distribution of protons and neutrons for 
^''O, ■*''Ca, ''*Ca and ^°^Ph are compared with the resuhs 
of N L3 (241] and experimental results of charge distributions 



the differences between the full HFB and the quasilocal 
approach plus BCS are less than 0.2%, which are similar 
to the values found for magic nuclei [68'|. 

In this paper we will use this formalism for the study of 
the binding energies of some open shell spherical nuclei. 
To this end we will use the mean field provided by the 
SEI. For a sake of simplicity we chose as pairing inter- 
action a density-dependent zero-range force of the type 
proposed by Bertsch and Esbensen [94| : 



1'(l"l,I"2) = Vq 



V 



Po 



5(ri-r2). (35) 



This effective pairing force is extensively used in nuclear 
structure calculations [l^, H^, [osj - looj . The numerical 
values of the parameters Vb=-481 MeV fm^, 77=0. 45, 
Q!=0.47 and po=0.16 fm~^ are taken from Ref. [sij . 
These parameters, together with a cutoff energy ec=60 
MeV from the bottom of the single-particle potential, 
were fitted in [84] to reproduce the gap values of the 
Gogny Dl force in neutron matter. We use here an 
improved BCS approach [23| where the resonant levels 
j92| are simulated very efhciently by quasibound single- 
particle energy levels retained by the centrifugal (neu- 
trons) or centrifugal plus Coulomb (protons) barriers. 
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FIG. 10: (Color Online) The two neutron (proton) and one 
neutron (proton) separation energies of Sn isotopes (N=82 
isotones) as a function of the neutron (proton) number in the 
uppe r (lower) panel in comparison with the experimental data 
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and charge radius Sr^h (lower panel) of 161 spherical even- 
even nuclei as function of mass number A. The experimental 
data for b inding energies and charge radii are taken from Refs. 

fioil.[Tol. 



The two-body center of mass correction has been taken 
into account self consistently by using a pocket for- 
mula based on the harmonic oscillator and derived in 
Ref. 100]. 
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FIG. 13: (Color Online) The average pairing gap from SEI for 
constant gap and constant strength approximation compared 
with the HFB prediction 



B. Results of finite nuclei 

In order to determine the three parameters, to, xq and 
Wo, we proceed as follows. Keeping the nine parameters 
determined in ANM (see Table I) fixed, (i) We reproduce 
the experimental binding energy of '^'^Ca by adjusting to, 
(ii) The spin-orbit strength parameter Wo is adjusted to 
reproduce the experimental splitting of the neutron and 
proton Ip levels in ^^O (iii) The parameter xo is obtained 
by fitting the binding energy and charge radius of ^'^^Pb. 
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FIG. 14: (Color Online) The odd-even mass staggering of 
Sn isotopes computed with the SEI energy density functional 
compares with the experimental values (upper panel). The 
same for isotones of N=82 (lower panel). 



The steps (i)-(iii) are repeated till the self-consistency of 
the parameters <oj and Wq is achieved. The value of 
the eleven parameters of the SEI entering in Eqn. ([1]) as 
well as the strength of the spin-orbit force are given in 
Table II. 

The binding energy per particle (BE/ A) and the pro- 
ton (rp), neutron (r„) and charge {rch — \J'''p^ ^.64 

fm) rms radii of the standard magic nuclei computed 
with the SEI quasilocal energy density functional are dis- 
played in Table III. These results are compared with the 
corresponding values obtained using the Gogny (DIS) 
[33„M] and the M3Y-P5 [M] forces and the RMF (NL3*) 
parametrization [2^ as well as with the experimental 
data (binding energies from [lOl| and charge radii from 
[l02]) that are also displayed in Table III. It is observed 
that the BE/ A obtained with the SEI functional matches 
nicely with the experimental data for these magic nuclei . 
In case of radii, our results are similar to those provided 
by the Gogny and the RMF calculations. The neutron 
skin thickness, defined as the difference between the neu- 
tron and proton rms, i.e. = r„ — r^, in ^"^Pb predicted 
by our calculation using the SEI functional is 0.21 fm. 
This value is in agree ment with the experimental result 
obtained in Ref. |103l | using proton-nucleus scattering. 

Experimental information about single-particle ener- 
gies of even-even nuclei can be obtained from the low- 
lying excited states of the adjacent odd nuclei by adding 
or picking up a single nucleon. The ordering of the par- 
ticle and hole levels, for both neutrons and protons, pre- 
dicted by t he SEI fu nctional is in agreement with the 
experiment (l04i Il05| except for the neutron 2/5/2 level 



which lies below the 3pi/2 one. This is, however, a rel- 
atively common fact in many mean field models, as for 
example the well reputed NL3 [106{ . The level spacing 
predicted by the SEI functional is larger than the ex- 
perimental one, but similar to the spacing obtained with 
the DIS and M3Y-P5 forces. The spin-orbit splittings of 
the single-particle levels 3p and 2/ of neutrons and 2d 
of protons obtained with our functional overestimate the 
experimental values as it also happens with the DIS and 
M3Y-P5 predictions. 

The radial dependence of the density distributions for 
protons and neutrons of the nuclei ^^O, '^'^Ca, ''^Ca and 
^°^Pb are displayed in the four panels of Fig.|9l The neu- 
tron and proton densities for the same nuclei obtained us- 
ing the RMF (NL3) |23| alongwith the exp erime ntal data 
for charge distributions, taken from Ref. [l07j |. are also 
given for comparison in the same figure. For ^^O, ''"Ca 
and ^^Ca the agreement between our densities and the 
ones predicted by NL3 model is quite good. The density 
distributions of ^°®Pb computed with SEI and NL3 dif- 
fer more between them, mainly in the bulk region. The 
proton densities predicted by the SEI functional agree 
reasonably well with the experimental charge distribu- 
tions at the surface. In the interior of the nucleus the 
quantal oscillations shown by our proton densities are, in 
general, well averaged. 

From Table III it is seen that the SEI functional pre- 
dicts ground-state energies and radii that accurately re- 
produce the experimental values. In this respect, the 
following comments are in order. To determine the nine 
parameters of ANM, we have considered the empirical 
values of three NM parameters, namely poj e(po) and 
Es{po). The widely accepted ranges of these parame- 
ters are po=0.17±0.03 /TO~^ e(po)=-16±0.02 MeV and 
Es{po)=30-35 MeV. We have examined the variations of 
these parameters within their ranges and have found that 
small variations in the values of pQ, e(po) and Es{po) 
have a large influence on the predictions of binding en- 
ergies and radii in finite nuclei, although the changes in 
the NM predictions are not of much significance. It is 
also found that po critically depends on the value of 7 
that determines the stiffness of EOS in SNM. For a given 
value of 7, there is a critical value of po for which the 
predictions of binding energies and radii will have min- 
imal deviation from experimental results. Thus for the 
EOS considered in this work corresponding to 7=1/2, the 
best results in finite nuclei are found for po =0-157 /m""^, 
e{po)= -16 MeV and Es{po) =35 MeV. 

We have performed further investigations to check the 
ability of our proposed SEI functional to describe ground- 
state properties of spherical nuclei. In this exploratory 
calculation pairing correlations are treated at BCS by 
the reasons pointed out before. Notice, however, that 
the BCS approach for pairing has been used in several 
well known mass tab les a s the ones computed with NL3 
[103 and HFBCS-1 [l0|. In the upper panel of Fig. M 
we display the two-neutron and the one-neutron separa- 
tion energies for Sn isotopes and in the lower panel of 
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the same figure the two-proton and one-proton separa- 
tion energies for the N=82 isotonic chain. To obtain the 
one-neutron and one-proton separation energies we have 
to compute odd-even or even odd-nuclei. To deal with 
these nuclei we perform, as in Ref. [s^ , spherical blocking 
calculations which neglect part of the core polaritazion 
effect in the odd nuclei fooj . In these calculations one has 
to block several single-particle states around the Fermi 
energy and search for the configuration giving the lowest 
energy. From Fig. [TUJ we see that the SEI predictions 
match fairly well with the experimental values for both 
chains. To be more quantitative, we have also computed 
the two- and one-nucleon separation energies in other iso- 
topic and isotonic chains with magic proton and neutron 
numbers, respectively. In Table IV we show the rms de- 
viation of these separation energies with respect to the 
experimental values |lOl| for the Z—2Q, 50, 82 isotopic 
and A''=50, 82, 126 isotonic chains. Our pred iction s are 
compared with the HFB results of Refs. jlll| and jll2l | 
(see the former reference for more details). A similar 
comparison between our results and the ones provided 
by the SLy4 and Gogny DIS forces and the NL3 param- 
eter set, but only for two-nucleon separation energies is 
given in Table V. From Table IV, it c an b e seen that our 
results are not so good than those of |lll| but similar to 
the o nes of [112]. Notice, however, that in Refs. |lll| and 
[ll2j the strength of the pairing force is fitted to some 
experimental data while in our case the pairing interac- 
tion is taken from the literature without any refit of their 
parameters to experimental data. The quality of the our 
predictions for two-nucleon separation energies also com- 
pare very satisfactorily with the results obtained from 
HFB calculations with SLy4 and SkM* [lli and DIS 



of charge radii, defined as Ar^^ 



' ch 



„2 (208 



[ll6j forces and the ones calculated with the NL3 pa- 



108 | 



rameter set 

We have performed more tests to check the ability of 
the SEI energy density for describing ground-state prop- 
erties of finite nuclei. To this end, we have computed the 
binding energy and charge radii of 161 even-even spheri- 
cal nuclei between ^^Ne and ^^^U. The difference between 
the theoretical prediction and the experimental value for 
binding energy /\E and charge radius Sr^h are shown in 
the upper and lower panel of Fig. 11, respectively. In 
general, the deviations in binding energies and charge 
radii for the 161-spherical nuclei are broadly reproduced 
within ±2 MeV and ±0.02 /to, respectively, baring few 
exceptions. The overall rms deviation in energy and 
charge radius of these nuclei are {/S.E)rms =1.5402 MeV 
and {l^R)rms = 0.0152 /m, respectively. These values 
are a little bit smaller than the corresponding rms devia- 
tions obtained with well calibrated effective interactions 
as DIS, SLy4, NL3 and BCP for the same set of nuclei 
(see Table 3 of Ref. [l^). We have also estimated 303 
binding energies and 111 charge radii of odd spherical 
nuclei, for which the experimental values are known. For 
these nuclei we find {AE)rms =1.6501 MeV and (Ai?)„ns 
= 0.0198 /to, respectively. 

Using our SEI, we have also computed the isotopic shift 



in Pb-isotopes. Analyzing the Ar^^ values obtained with 
the SEI, one can see in Fig. [T^ that they lie quite close to 
the experimental values [102] and that the kink at ^=208 
is reasonably well reproduced. It is known that Skyrme 
and Gogny forces use isospin independent spin-orbit in- 
teractions those are unable to reproduce the experimen- 
tal kink exhibit by the isotopic shift of the charge ra - 
dius in Pb isotopes at the double magic 208pb pl3iill4j. 
This mismatch can be cured by introducing an isovec- 
tor spin-orbit contributi on as it is the case of the Ski 
family of Skyrme forc es 114l | or using non- linear RMF 
parametrizations |ll7j . However, the SEI functional in- 
corporates a spin-orbit contribution that does not con- 
tain isovector part (see Eqn. (j23p ) but clearly shows up 
the kink without imposing it as a constraint in the fit- 
ting procedure of the functional. A possi ble explantion 
is the following. In a recent paper jll8| it is claimed 
that the development of the kink in the isotopic shift of 
the charge radius of ^°®Pb is largely determined by the 
occupation of the liii/2 neutron level. This orbital has 
a principal quantum number n=l and overlaps strongly 
with the majority of the proton orbitals. This produces a 
relatively large pulling of the proton orbitals by the neu- 
tron ones via the symmetry energy increasing the proton 
radius, and, therefore, developing the kink at A — 208. 
From Fig. 8 we can see that the SEI functional predicts 
that the ljii/2 lies close to the 2gg/2 one and the same 
happens for heavier Pb isotopes. Therefore, one can ex- 
pect a relatively large occupation of liii/2 state. This 
fact together with the relatively high symmetry energy 
of the SEI seems to be the reason for the appearence of 
the kink exhibit by the SEI results in Fig. 12. 

We shall now discuss in some detail the pairing proper- 
ties of our SEI energy density functional. This functional 
includes a pairing contribution coming from a zero-range 
density-dependent force that simulates the Gogny inter- 
action in the particle-particle channel. As far as the 
Gogny pairing force gives a good description of finite 
nuclei combined with a reasonable mean field, not only 
provided by the Gogny interaction but also by Skyrme 
forces and RMF parametrizations [1^, we are, actu- 
ally, checking the single-particle energies obtained with 
the SEI energy density functional. In Fig. [13] we dis- 
play the average gaps along the whole Sn isotopic chain 
obtained with our formalism compared with the values 
predicted by a full HFB calculation with the Gogny 
DIS force in both, particle-hole and particle-particle, 
channels. From this figure it can be seen that the av- 
ejaged gaps, defined as A^,2 = J2n'"l^ri/ J2n^l ^^^^ 
A„i, = X]„ UnVn^n/ J2n "n^" uiccly rcproduce the HFB 
values in this isotopic chain taken of Ref. [o^l . This result 
suggests that our model, with a mean field part whose 
parameters are fitted to ANM and magic nuclei together 
with a realistic pairing force taken from the literature, is, 
in principle, well suited for describing open-shell nuclei at 
mean field level with a quality similar to that found us- 
ing well known effective interactions or RMF paramtriza- 
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tions. 

A possible way of estimating pairing correlations from 
an experimental point of view is throu gh the so-called 
odd-even mass staggering (see e.g. Ref. f90'|). In Fig. [HI 
we display the odd-even mass staggering using a three- 
point formula [o^] 



A. 



E{N + 1) - 2E{N) - E{N ~ 1) 



(36) 



for the Sn isotopic chain (upper panel) and for the iV=82 
isotonic chain computed with the SEI energy functional 
in comparison with the experimental values. Our the- 
oretical prediction of the mass staggering slightly over- 
estimate the experiemental values in both, Sn isotopic 
chaim and A^=82 isotonic chain. Therefore, some com- 
ments are in order. On the one hand, the odd-even mass 
staggering of experimental masses is not a pure measure 
of the pairing correlations [l^l ■ It also contains mean field 
contributions related to the rotational and time-reversal 
symmetries breaking not accounted by the simple block- 
ing approach. On the other hand, the fact that the theo- 
retical odd-even mass staggering overestimates the corre- 
sponding experimental values may be not so dramatic. It 
is known that when pairing is computed using the Gogny 
interaction, the average gaps are larger than the mass gap 
[33 |. However, as discussed in this reference, effects be- 
yond mean-field, such as quasiparticle vibration coupling, 
are expected to reduce the average gaps. Finally to point 
out that in this exploratory calculation no attempt has 
been done in order to adjust simultaneously pairing and 
mean field. Therefore, it seems reasonably that a slightly 
different protocol in the fitting procedure of the param- 
eters to and xq, for example including some open shell 
in it, would allow to find an odd-even mass staggering 
in better agreement with the experimental data. From 
all the discussion developed along this section, we believe 
that the reported calculations show up clearly that the 
proposed energy density functional based on the SEI is a 
realiable tool for dealing, at least in spherical nuclei, with 
quality similar to that found using well calibrated Gogny 
or Skyrme forces or succesful RMF paraametrizations in 
the same scenario. 



V. SUMMARY AND CONCLUSIONS 

We have proposed a simple effective interaction aimed 
to describe accurately the main trends of microscopic cal- 
culations in nuclear and neutron matter and, at the same 
time, to reproduce the ground state properties of finite 
nuclei with a quality similar to the one obtained using 
successful effective interactions of Skyrme, Gogny and 
M3Y type or relativistic mean field model. Our interac- 
tion contains a single Gaussian form factor for describ- 
ing the finite range part of the force plus a zero-range 
part that includes a density dependent term to simu- 
late the effect due to three-body forces. Our effective 
interaction depends on eleven parameters, nine of them 



fitted to nuclear and neutron matter data. The SNM 
is fully determined by six parameters, which are three 
strength combinations, the range of the force, the expo- 
nent 7 and the factor b. Two of these six parameters, the 
strength of the exchange interaction and the range of the 
force, completely determine the momentum dependence 
of the mean field in SNM. These two parameters are con- 
strained by using the optical potential data which im- 
poses that the calculated mean field in normal SNM van- 
ishes at a kinetic energy 300 MeV of the incident nucleon. 
The exponent 7 of the density dependent term is chosen 
amongst the values that gives the pressure-density rela- 
tion in agreement with the experimental results of HIC 
and production. The parameter b is fixed to avoid 
the supr aluminous behaviour in NM. The two remaining 
strength combinations in SNM are determined from the 
saturation conditions. In our model, the range of the 
force is chosen to be the same for interactions between 
pairs of like and unlike nucleons. Under this restriction, 
the splitting of the exchange strength parameter in ANM 
is decided from the condition that the entropy in PNM 
does not exceed the one in SNM. The splitting of the two 
other strengths into like and unlike channels in ANM are 
decided from the value of the symmetry energy and by 
imposing that the nucleonic part of the energy density in 
charge neutral /3-stable matter be a maximum. With the 
nine parameters of our effective interaction determined in 
this way, the corresponding EOS and mean field follow, 
quite closely, the general trends shown up by sophisti- 
cated microscopic calculations. 

To describe finite nuclei we use this interaction to build 
up an energy density functional in the framework of the 
quasilocal Density Functional Theory. In particular we 
use the semiclassical h?- expansion of the density ma- 
trix to deal with the exchange contribution. This proce- 
dure allows to write the single particle equations in finite 
nuclei in a similar form as in the case of the Skyrme- 
Hartree-Fock equations. We determine the open param- 
eters of our interaction to and xq , and the strength of the 
spin-obit force Wq by a simple fit to experimental data 
of three closed sheU nuclei ^^O, "OCa and ^ospb. With 
this choice the binding energies and charge radii of the 
standard magic nuclei are accurately reproduced. With 
our simplified interaction fully determined, we explore 
its predictive power for describing open-shell spherical 
nuclei. To this end, we add to our functional the con- 
tribution of a zero-range pairing interaction, taken from 
the literature, which reproduces the Gogny neutron gaps 
in PNM. In this exploratory calculation we use the BCS 
approximation instead of HFB theory for dealing with 
open shell nuclei and the simplified spherical blocking 
approach to estimate the ground-state properties of odd 
nuclei. It is important to note that within this frame- 
work, the results for open-shell nuclei are predictions of 
our model as far as no parameter has been fitted to open- 
shell nuclei. We have computed the energy of 161 and 
charge radii of 88 spherical even-even nuclei in the mass 
region A = 16-224. Our calculation reproduce the ex- 
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perimental values fairly well. The rms deviations of the 
binding energies i and radii of these nuclei with respect 
to the experimental values are 1.5402 MeV and 0.0152 
fm, respectively, which are slightly better than the same 
results for the set of nuclei computed using standard in- 
teractions such as Gogny DIS, BCP, NL3 and SLy4. We 
have also calculated the one- and two-nucleon separation 
energies along the Ca, Sn and Pb isotopic chains and 
7V=50, 82 and 126 isotonic chains. With our model, the 
rms deviation of the theoretical predictions respect to 
the experimental data compares favourably with the re- 
sults obtained using well calibrated effective interactions 
of Skyrme, Gogny and relativistic mean field type. We 
have also investigated the isotopic shift of the charge- 
radii in Pb. We find, surprisingly, that the kink shown 
by the experimental charge radii at ^=208, is qualita- 
tively reproduced by our model, in spite that our spin- 
orbit force has no isospin dependence, which is needed, 
in principle, to reproduce the kink with non-relativistic 
mean field models. We have given a possible explanation 
of this fact based on the small gap between the 2gQ /2 and 
liii/2 neutron levels and the relatively large symmetry 
energy predicted by the SEI model. 

In this exploratory analysis of finite nuclei properties 
described with our simplified effective interaction we have 
restricted to spherical nuclei. To improve the treatment 
of pairing correlations by replacing the BCS approach 
by a Hartree-Fock-Bogoliubov calculation, and extend 



our study to deformed nuclei and some collective excited 
states is a necessary task to confirm the success of our 
simplified effective interaction that will be developed in 
future works. 
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